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Abstract. 

Wilson loop averages of pure gauge QCD at large N on a sphere are 
calculated by means of Makeenko-Migdal loop equation. 
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The exact solvability of QCD2 was noticed long time ago by A.A.Migdal 
0. To realize it, one can integrate over U{N) link variables using only the 

orthogonality of group characters. 
Following this, an exact solution of QCD2 on an arbitrary 2-manifold was 

obtained in ref. [@]. 

In particular, the partition function for a closed surface of a genus g and of 
area A (coupling constant absorbed into area) has the form: 

ZM) = T.df-'^^w{-^Cn) (1) 

where dji is the dimension of an R-th irreducible representation , 

n (1+^) (2) 

l<i<j<Af J ^ 

Cfi is the value of the second Casimir operator, 

N 

CR = J2niin, + N + l-2t) (3) 

where rii are highest weight components of i? = {rii > n2 > . . . > n^v}. 

The Wilson loop average has the form P]: 

M A 

w,{c)= E *«.-Mn4r^^^-p(-|^c,,) (4) 

Ri,---,Rm k=l 

where M is the number of windows, A^s are their areas, Qk is the "genus 
per window" and the coefficient ^r^...rm is the U{N) group factor related to 

a form of the given contour Q (see ref. for details). 

Recently, the partition function (|I|) has been reexpanded over and at 
large N it has been interpreted as a sum over branched coverings 0. 
In order to develop a large N technique appropriate for calculations, say, of 
loop averages, it is tempting to rewrite the partition function (|ID at large N 

^ In particular, for a contour without self- intersections , this factor is the product, taken 
over all pairs of neighboring windows < ij >, of Wigner coefficients ^r-r^ (see eq.(p^ 
below) . 

^The similar interpretation has been known for the Wilson loop operators as well B. 
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as a path integral over continuous Young tables ^ . In the paper of Douglas 
and Kazakov 0, in the framework of this approach, a third order phase 
transition on a sphere was found. Recently, Boulatov has calculated the 
Wilson loop average for the simplest contour [0]. The goal of present paper 
is the calculation of an arbitrary Wilson loop average at large N on a sphere. 
First, we repeat briefly the steps of In the spherical case, g = 0, the 

sum becomes divergent at small areas. Hence, in the large limit, a 
non-trivial saddle-point should exist. Let us introduce the continuous 

function 

h{x)=lun^^{^-^-n,); ^ = ^'1 (5) 
Then the saddle-point equation takes the form 
\ = (6) 

Here, 



2 J h-y 



obeys the inequality 
p{h) < 1 (8) 
If eq.(|]) is ignored, then the solution of eq.(|^) is the semi-circle distribution, 

(9) 

which is valid for the areas A < tt"^. For the areas A > tt"^, the inequality (|^) 

is crucial and the solution of eq.(|^) is 




-^H^i^'-h'W-mi{-^,^) for -a <h<-b 
p{h) = { 1 ioi -h<h<h (10) 



^ J(a2 - h-W - &^)ni for b<h<a 
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where Ili{x, k) is the complete elhptic integral of the third kind with the 
modulus k = ^ and parameters are to be determined from the equations 



a{2E - k'^K) = 1 aA = AK (11) 

At the critical value, = tt^, the third order phase transition takes place. 

Since the saddle-point solution is known in both phases, one can calculate 

the Wilson loop averages (^) at large N on a sphere. 
The simplest Wilson loop on a sphere corresponds to M = 2 and 

91 = 92 = ^ (see (§), 

where 

^ns = J dUXn{U)Xs{U^)tT{U) (13) 

is the multiplicity of an irreducible representation S in the tensor product 
of R with the fundamental representation f; A = Ai + A2 is the total area. 

The expression for this quantity valid for both phases is [0]: 

WiA,,A2) = i — e^^'^-^C^) = / ^e-^^'^+^C^) (14) 
Jc 2'Ki Jc 27ii 

where 

and the contour C encircles the cut of f{h) (contour C encircles the cut in 
the opposite direction). Obviously, W{Ai,0) = W{0,A2) = 1. 

In the weak coupling phase, f{h) = — \J —~A, and in terms of the 

new variable z = ^ the integral takes the form 

,y„.(A.A.)^-i=£i|(l.l)e— I (16) 
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where a = ^ and (3 = Expanding the exponential and taking the 

residue at zero, one find 



„ oo 
•^<^ n=0 



(17) 



(-1)" (n'f^\n _ / A1+A2 J ( /TaT^ 



1A2 
-A2 

n=0 



In what follows we will also need A ^ 00 limit of (14). As it has been 
shown in ref. 0, the function (^) in this case takes the form 
/(x) = log If^Y- -^^^ example, in the case = 00, we have 

2 Jo Zm X + 1 

i.e., standard QCD2 area law. 
In principle, all set of (P can be calculated in the same fashion. However, 
since an expression for simple loop ([T^) is known, all others Wilson loops 

can be derived easily from the loop equation. 
The equation for Wilson loops derived by Yu.Makeenko and A.Migdal a 
long time ago P, has especially simple form in two dimensions [^j: 

RiWn{C) = W,{Cl)Wn-,-i{Cl) I = l,...,n (19) 

where n is a number of self-intersections of contour C. Each equation of 
system (|^) corresponds to given i-th. point of self-intersection of C {C} and 
are results of disconnection of contour C in this point)^. Here, 
differential operator Ri acting in i-th self-intersection point is the linear 
combination of derivatives with respect to areas of windows ii, 
neighboring at this point (in what follows, we use the notation dk = j^). 



Ri = di^ - (9i2 + - di^ (20) 

where different signs correspond to windows having common boundary ^ 

(Figi). 



^ Actually, on a sphere (as well as on any compact surface), the number of independent 
equations of the system ( |l9|) can be less, than n. Nevertheless, it turns out that in the 
spherical case the system (|19|) is sufficient to derive a complete set of loop averages. 

''We do not fix the total sign here and below since it is always can be restored from the 
condition VK(0,...,0) = 1. 
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Figure 1: Disconnection of a contour at the self-intersection point 



If all Wilson averages for loops with number of self-intersection less, than n, 
are known, then we know an explicit form of RHS of (|19]). Solving 
recurrently the system (|^) we get all set of loop averages. 
Starting point of (HJ) is equation for 8-Iike contour drawn on Fig.2. In 
this case, the RHS of (|1^) has the form of product of two simple loops, Wq 
(eq.(|T^)), which could not be derived from the Makeenko-Migdal equation. 

The equation, corresponding to Fig.2, has the form 



RW,{A^,A2,A,) =f — f -^e^'^^+'^'y-^^^^-^^y) (21) 
J Zm J Zm 

where 

R = di + d2- 2^3 (22) 

The natural anzatz is: 

WMi, A,)=l^l p_,,Ar.+A,y-Kx)-mF (23) 
J Zm J Zm 

Then, function F = F{x, y; Ai, A2, A3) should satisfy the equation: 
{R + x + y)F =1 . (24) 

The solution of (0) is: 

^ ^ l + e-(-+^)('-^+^)e(Ai,A2,A3) 

x + y ^ ^ 

where function ^{Ai, A2, A3) should satisfy the equation: 
= (26) 
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and c is arbitrary constant. An explicit form of function ^ and value of 
constant c should be defined from obvious conditions: 

W^l(0,/l2,A3) =W^o(^,^3) , (27) 

0,^3) =1^0(^1,^3) , (28) 
which fix unambiguously the function F. The result is: 

WMuA..A,).fl^J^^—^— (29) 

Let us notice that this expression can be equally rewritten as 
WMuA.,A,)^pf-J^; — (30) 

where P means that integration over x should be performed after 
integration over y. The equivalence between (^) and ( PP] ) corresponds to 
the fact that the contour drawn on Fig. 2 coincides with one drawn on Fig. 3. 
The known answers for a plane |^-[|ri|] are follow from ( PU] ) and ( pOD as well: 



1^1(^1,^2,00) = e-^^ (31) 
1^1(^1,00,^3) =e-^^-^(l-Ai) (32) 



A 



Wi(oo,A2,A3)=e-^^--(l-A2) . (33) 

In the case of two self-intersections there are two different types of 

contours (Fig.4). 

Since Wq and Wi are known, we are able to write down explicitly an 
equation for W2^s. For example, for contour Fig.4 (a), we have 

i?i 1^2(^1, ^2, ^3, A,) = WoiA,, + + ^4)1^1(^2, A3, Ai + A,) (34) 

R2W2{Au A2, A3, A4) = 1^0(^3, A^ + A2 + ^4)1^1(^1, A2, A3 + A4) (35) 

where 

^1 = ^1 + ^2 - 2^4 , (36) 
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R2 = d2 + d;-2di . (37) 
Repeating the steps of previous case, we derive 

while for contour Fig. 4 (b), the answer is 
f^e-/(-^)-ff^e-/(-3) ^, ^ .. ^ (39) 

It is also instructive to consider one more example, the contours where pair 
(or pairs) of windows are separated by more than one self-intersection 
point. The first example of such a contour we match while consider the case 
of three self-intersections. The contour drawn on the Fig. 5 has two 
windows, 1 and 4, separated by three self-intersection points. In this case, 

the solution of loop equation is 

Ws{A,,...,A,) = 

(40) 

r d2L r dy ^- f(x)~ f(v) e^^-+Mv / A.fa+ty) , ^Asix+y) _ ^(A2+A3)(x+y)\ 
•f 2iTi * 2iTi x+y ^ ' 

Let US notice that the structure of expression under integral in (|iO| ) is in 
one to one correspondence to structure of the plane solution which we 

easily reproduce from 



Vr3(^i,-,^4,oo) =e-^(e-^2 +e-^3 _g-A.-A3) ^4^) 



Now, following previous consideration, we are able to formulate general 

solution for arbitrary Wn[C\. 

i e^.-.-/(x.) -Q J2 (-l)5e-^"("»±"^) 

27rz J 2m j^^^ ^-^^ Xi ± xj ^^-^ 

(42) 

where: m is number of windows of the disk topology appearing as a result 
of disconnection of C at all self-intersection points (we associate variables 
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Xk with all such a windows) Q Ak is an area of such a windowQ the second 
product goes over all pairs < ij > of windows from this set which are 
separated by (at least one) self- intersection point; minus sign of 
"propagator" ^ corresponds to the case of self-overlapping, when i-th 
window is internal with respect to j-th one; P-symbol means that, in the 
latter case, integration over xj should be performed before integration over 
Xi] if i-th and j-th windows are separated by riij > 1 self- intersection points 
Q, then the sum goes over all disconnections at such points and A]^ is the 
area of the overlapping region appeared as result of g-th disconnection (see 

eq.(ig)). 

Eq. (^) is the answer for arbitrary Wilson loop average in large N pure 
QCD on a sphere. Let us emphasize that it is valid in both weak and 
strong coupling phases. A complete information about phase transition is 
contained in the function p{x) {f{x) as well). This function is given by 
eq.(^) in the weak coupling phase and by eq.(|TDp in the strong coupling 

phase. 
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